% 
chang = 3.41-2* 0.275;
banjing = 16*0.55;
jiaodu = 1/banjing;
x_tou = zeros(2,101);
qiyu_zuobiao1 = zeros(224,101);
qiyu_zuobiao2 = zeros(224,101);
jiaodu_quanbu = zeros(101,1);
for i = 350:450
    jiaodu1_huchang = 16*2*pi-jiaodu*i;
    jiaodu_quanbu(i-349,1) = jiaodu1_huchang;
end
%计算头的位置坐标
for i =-100:0
    jiaodu1_huchang = jiaodu_quanbu(-(i-1),1);
    x_tou(1,-(i-1)) = 3*cos(jiaodu1_huchang);
    x_tou(2,-(i-1)) = 3*sin(jiaodu1_huchang)+3;
    leiji_theta = jiaodu1_huchang;
    bianjingwenti = 3/(jiaodu1_huchang)*jiaodu1_huchang/(2*pi);
  for j = 2:223   
      bianhua_heta = 1.65/bianjingwenti;
      leiji_theta = leiji_theta+bianhua_heta;
      bianjingwenti =  3/(jiaodu1_huchang)*leiji_theta/(2*pi);
     qiyu_zuobiao1(j,-(i-1)) = bianjingwenti * cos(leiji_theta);
     qiyu_zuobiao2(j,-(i-1)) = bianjingwenti * sin(leiji_theta);
  end
end
x_tou = zeros(2,30);
qiyu_zuobiao1 = zeros(224,30);
qiyu_zuobiao2 = zeros(224,30);
jiaodu_quanbu = zeros(30,1);
for i = 525:550
    jiaodu1_huchang = 16*2*pi-jiaodu*i;
    jiaodu_quanbu(i-524,1) = jiaodu1_huchang;
end
%计算头的位置坐标
for i =0:30
    jiaodu1_huchang = jiaodu_quanbu(i+1,1);
    x_tou(1,i+1) = 3*cos(jiaodu1_huchang);
    x_tou(2,i+1) = 3*sin(jiaodu1_huchang)+3;
    leiji_theta = jiaodu1_huchang;
    bianjingwenti = 3/(jiaodu1_huchang)*jiaodu1_huchang/(2*pi);
  for j = 2:223   
      bianhua_heta = 1.65/bianjingwenti;
      leiji_theta = leiji_theta+bianhua_heta;
      bianjingwenti =  3/(jiaodu1_huchang)*leiji_theta/(2*pi);
     qiyu_zuobiao1(j,i+1) = bianjingwenti * cos(leiji_theta);
     qiyu_zuobiao2(j,i+1) = bianjingwenti * sin(leiji_theta);
  end
end
qiyu_zuobiao1 = zeros(224,75);
qiyu_zuobiao2 = zeros(224,75);
jiaodu_quanbu = zeros(75,1);
for i = 450:525
    jiaodu1_huchang = 16*2*pi-jiaodu*i;
    jiaodu_quanbu(i-449,1) = jiaodu1_huchang;
end
%计算头的位置坐标
for i = 0:75
    jiaodu1_huchang = jiaodu_quanbu(i+1,1);
    x_tou(1,i+1) = 3*cos(jiaodu1_huchang);
    x_tou(2,i+1) = 3*sin(jiaodu1_huchang)+3;
    leiji_theta = jiaodu1_huchang;
    bianjingwenti = 3/(jiaodu1_huchang)*jiaodu1_huchang/(2*pi);
  for j = 2:223   
      bianhua_heta = 1.65/bianjingwenti;
      leiji_theta = leiji_theta+bianhua_heta;
      bianjingwenti =  3/(jiaodu1_huchang)*leiji_theta/(2*pi);
     qiyu_zuobiao1(j,i+1) = bianjingwenti * cos(leiji_theta);
     qiyu_zuobiao2(j,i+1) = bianjingwenti * sin(leiji_theta);
  end
end